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Abstract. We study the approximation properties of a harmonic function 
u £ k > 0, on relatively compact sub-domain A of Q, using the Gen- 

eralized Finite Element Method. For smooth, bounded domains Q, we obtain 
that the GFEM-approximation us satisfies ||u — u s\\h 1 (A) ^ CYO IMIffi-fc(o.) > 
where h is the typical size of the "elements" defining the GFEM-space S and 
7 > is such that the local approximation spaces contain all polynomials of 
degree k + 7 + 1. The main technical result is an extension of the classical 
super-approximation results of Nitsche and Schatz [?] and, especially, [?]. It 
turns out that, in addition to the usual "energy" Sobolev spaces H 1 , one must 
use also the negative order Sobolev spaces H , I > 0, which are defined by 
duality and contain the distributional boundary data. 
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Introduction 

To motivate our results, let us consider the boundary value problem 

J Am = on Q, 

\ % = g on T := dtl, 

where is a smooth, bounded domain in R n , <9fi is the boundary of f2 and g £ 
H r -*/ 2 {d£l), r g E. The case r > 3/2 was extensively studied. Here we are 
interested mainly in the case r < 3/2. We are looking for a solution u g H r (Q). 

For r > 3/2, the boundary values (or traces) u\qq and d„M|on are defined clas- 
sically, because the restriction to the boundary extends by continuity to a map 
H r (Cl) 3 u -> d v u € H r ~ 3 / 2 (dn), see [?, ?] for example. For r < 3/2, this is no 
longer true, but then one takes advantage of the fact that u satisfies an elliptic 
equation, so it is still possible to define d v u g H r ~ 3 / 2 (dQ), see [?, ?]. We can 

I. Babuska was partially supported by NSF Grant DMS 0341982. V. Nistor was par- 
tially supported by NSF Grants DMS 991981 and 0200808. Manuscripts available from 
http://www.math.psu.edu/nistor/. 
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assume, without loss of generality, that f2 is connected. It is not difficult to prove 
that a solution u of Equation (JJ exists for any g such that (g, 1} = and that this 
solution satisfies 

( 2 ) IMIff''(tt) < C\\9 llffr-3/ 2(an) , 

with a constant C that may depend on r but is independent of g. This result will 
be discussed in detail in [?], where more references will be given. (This result will 
be used in this paper only as a motivation for our work.) 

Recall that A <g f2 means that A C £1 is bounded and dA and 90 are disjoint. 
If A <s is an open subset, then u satisfying will be smooth on A, regardless 
of what r is, and 

(3) \\u\\ H ^{A) <0\\u\\ H r( d Q), 

with a constant C that depends on A, Q, r, and m, but is independent of u satis- 
fying Aw = 0. An important problem, with potential practical applications, is to 
approximate on A the solution u of Equation Q . 

In this paper, we prove several results on the approximation of the solution u on 
subsets A <s f2. The main result of the paper, Theorem 14. 81 then gives in particular 
that the solution u of the boundary value problem Q satisfies 

(4) ||u - u s \\ H i(A) < C7i 7 ||u||iri-k(n), 

where us £ S is the Galerkin approximation and S is the Generalized Finite Ele- 
ment Space associated to a partition of unity {4>j} subordinated to a covering {u>j} 
of f2 satisfying Assumptions A-D of Section |2 provided that our local approxima- 
tion spaces contain all of polynomials of degree 1 + 7 + fc, 7 > 0. We stress that 
our results require not just the energy Sobolev space H 1 , but also negative order 
Sobolev spaces H~ l , defined in this paper as the duals of H, I G Z+. One of the 
main reasons for the need to considere the negative order Sobolev spaces is that 
the solution u is not in but in H 1 ~ k (U), in general. Moreover, even if we 

approximate the boundary data g and the solution u with functions in H , then 
the norm on a negative order Sobolev space will still have to appear in the estimate 
of the error. 

Here is now a brief description of the contents of the paper. We continue to 
assume that Q is bounded and connected, but we no longer assume that ft is smooth, 
except when explicitly mentioned. In Section ^ we set up the notation and we 
establish our conventions on Sobolev spaces. Section contains a quick review of 
the necessary definitions involving the Generalized Finite Element Method (GFEM) 
and their variants used in this paper [?, ?, ?]. This section also contains the 
assumptions that we make on our covering and partition of unity used to define the 
GFEM-space S. The space S will contain our approximate solution to the boundary 
value problem. The following section, Sectional contains the calculations necessary 
to establish our interior estimates. Our approach follows, to a certain extend, that 
in the article of Nitsche and Schatz [?] , relying also from Wahlbin's survey article 
[?]. The main differences between our paper and [?, ?] are due mostly to the fact 
that several assumptions from these papers are not fully satisfied in our approach. 
As in these articles, the main step is a super-approximation property, Proposition 
13.61 The proof in [?, ?] cannot be used to obtain Proposition 13.61 because the 
property "d a u = if \a\ is large," is not satisfied in general for u E S. For the 
results of Sectional we assume that f2 is smooth. Then we extend the definitions of 
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the Galerkin approximation us E S and of the form B(w, v) :— J n Vtifi) • Vv(x)dx 
to the case when v E H k ^ 1 (^l) is arbitrary and u E H 1 ~ k (il) can be written as 
u = M1 + U2, where Aui = in distributions sense and 1*2 E if (0). In what follows, 
1 — fe will play the role of r above. Several estimates for u and its approximation 
us are established in this section, including the main theorem, Theorem l4. 81 (whose 
main conclusion was repeated in Equation Q above). The last section, Section [5] 
contains a proof that, for a domain with piecewise C 1 -boundary, we can construct 
a family of partitions of unity with typical size of supports h — ► that satisfies 
Assumptions A-D for a fixed choice of the other values of the parameters (i.e., of A, 
B, Cj, k, m, and a). For this construction we assume that the local approximation 
spaces are *f?j = Qk, the space of polynomials of degree at most k. In particular, 
the various assumptions made in the results proved in the preceding sections, are 
satisfied for this family of partition of unity, and hence our results are not empty 
for domains with piecewise smooth boundary. By contrast, it is not possible to find 
a family of partitions of unity as above for domains with cusps, see Remark 15.71 
For suitable g, we plan to perform some concrete numerical simulations in a future 
paper. 

We shall write x := y if x is defined by y. 

The second named author thanks G. Grubb, A. Schatz, and L. Wahlbin for some 
useful references. A. Schatz has also made some useful comments on an earlier 
version of the manuscript, for which we are greatful. 

1. Preliminaries 
We begin by fixing the notation and terminology. 

1.1. Preliminary notation. We denote by R the set of real numbers and by 
C := {a + bi,a,b E R} the set of complex numbers. Also, N = {1,2,...} and 
Z + = {0} U N. Let x ■ y := x±y\ + X22/2 + ■ • • + %nVn be the inner product of two 
vectors x,y E M™. We shall denote by 

(5) f(y) := / e^y x f(x)dx 

the Fourier transform of /, as usual. By L 2 {VL) we shall denote the space of square 
integrable functions / : f2 — > C, for any domain £1 C R™. 

Let s > 0. Then iP(R n ) is the space of functions / £ L 2 (R") such that 

(6) II/IIW") :=(2^)-"£(i + |y| 2 n/(y)| 2 %<^. 

Here |y| = ^yf + .-. + yl if y = (yi, . . . ,y n ) E R". Let Q c R™, then H S (Q), 
s > 0, denotes the restrictions to 17 of functions / E H s (W l ), that is, 

(7) H S (Q) := {f\n, I E i7 s (R")}. 

The norm on H s (fl) is then ||/i||ij«m) := inf ||/||jy a (K")) the infimum being taken 
over all functions u E H s (M. n ) such that f\n = h. When s is a positive integer and 
f2 is a nice domain (Lipschitz, for example), the norm ||w||//s(f2) is equivalent to the 
usual norm 

(8) \\h\\ 2 s = II^HW), 

|q I <s 
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where a = (a\, . . . , a n ) G Z™ , \a\ := ce\ + 012 + . . . + a n , and d a :~ d^d^ 2 ■ ■ ■ d" n , 
as usual. The space Hq (f2) is defined as the closure of (ft) in Hq (f2) . 

1.2. Distributions. Let B R (0) denote the open ball of radius R centered at the 
origin. Also, let C£°(R™) be the set of infinitely differentiable, complex valued 
functions that vanish outside the ball Br(0), for some large R > 0. A linear map 
u : C C °°(R") -> C is called a distribution on R n [?, ?, ?] if, for any R > 0, there 
exists m G Z + = {0, 1,2,.. .} and C > such that 

(9) |«(0)| <C ^ \\d a 4>\\ L ™, if </> G C C °°(R") and = outside B fl (0). 

a I <m 

This definition does not exclude the case when larger and larger values of m have to 
be chosen as R — + 00, and in fact this situation actually occurs in specific examples. 
The set of distributions on R n will be denoted V'(R n ). 

We now fix more notation and terminology. If / is a function, then the closure of 
the set {/ ^ 0} is called the support of f and will be denoted supp(/). We shall also 
write (u, <fr) := u{4>) for the value of the distribution u on the function <f> G C£°(R"). 
Therefore, any <j) G C^°(U) has compact support. The support of a distribution u 
is the smallest closed set F such that (it, <j>) = for any (f> G C^°(R n \ F). 

Here are some examples of distributions. 

Example 1.1. If / is a measurable function on R™ that is integrable on any closed 
ball in R™ (i.e., it is locally integrable, or / G L 1 1 oc (R™)), then we can define 

(10) (/,<£) := / f(x)<t>(x)dx, 

for any G C C °°(R"). Thus any / G LlJW 1 ) defines a distribution on R™, that 
is, L 1 1 oc (R") C P'(R"). In this situation we shall say that our distribution / is, in 
fact, a locally integrable function, or, that our distribution is defined by a locally 
integrable function. 

The following two examples are relevant for the discussion of concentrated loads 
and moments. 

Example 1.2. The Dirac measure at a G R™ is the distribution 5 a defined by 

(5 a ,4>) := 4>(a). 

An explicit calculation shows that S a G i/~ n / 2 ~ e (R") and ||i5 a || ff -r 1 /2- £ ( R r 1 ) — > 00 as 

Example 1.3. The derivatives d a u of a distribution u arc defined by 

(d a u,(t)) := (-l)\ a \u,d a (t)). 

We now define the negative index Sobolev spaces. Thus, the space H~ s (R n ), s > 
0, consists of all the distributions u G X>'(R") such that there exists h G Ll oc (M. n ) 
satisfying 

(«,/)=/ h(y)f(y)dy, for all / G C C °°(R") 

and 

(11) := (27T)"" / (1 + |y| 2 ) s |My)| 2 ^ < 00. 
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The reader will recognize that the condition above is analogous to the condition of 
Equation ©. The difference is that now we allow s to take on negative values as 
well. 

The following alternative definition of the negative order Sobolev spaces will 
be useful later on. One first checks directly using the Fourier inversion formula 
together with Plancherel's formula that 

u(x)<f)(x)dx 



(12) \\u\\± s :=mf 



where u, <jy £ C c °°(R n ), v ± 0. 



Then H~ s (W l ) can be defined as the completion of C£°(R") in the norm given by 
the formula of Equation l|12l) . In particular, we obtain that H~ s (M. n ) is canonically 
isomorphic to the dual of H s (W l ). 

Similarly, for any open subset H C R™, we define the space H~ s (Sl), s > 0, as 
the dual of H S {Q). 

HIh-* ( o)=su P jfr'^ 1 , 0^etf s (fi). 

11011^(0) 

Our definition of negative order Sobolev spaces by duality follows [?, ?, ?], for 
example. Note however, that the negative order Sobolev spaces are often also 
defined by restriction from R™, as in [?, ?, ?], for example. The space of restrictions 
to Q of distributions in H~ s (R n ) is the dual of Hq(Q), and will be denoted Hq S (Q). 
The spaces Hq s (CI) = Hq(SY)*, s > 0, will also be used below. 

When il = R™, these two approaches yield the same spaces, but for general 
they may lead to different "negative order" Sobolev spaces. 

2. The Generalized Finite Element Method 

We now recall a few basic facts about the Generalized Finite Element Method 
[?, ?, ?]. This method is especially convenient since it provides finite element spaces 
with high regularity. Most of the results of this section work for a general bounded 
open set fi, except the application to the boundary value problems, Subsection 12. 21 

2.1. Basic facts. We shall need the following slight generalization of a definition 
from [?,?]: 

Definition 2.1. Let SI C M™ be an open set and {uij} be an open cover of 
with no k + 1 of the sets u>j having a non-empty intersection. Also, let {<fij} be a 
partition of unity consisting of Lipschitz functions and subordinated to the covering 
{u>j} (i.e., suppt^j C ujj). If 

(13) \\d a 4> 3 \\ L ^ m <C k /(dmmuo j )\ k = \a\<m, 

then {ujj} is called a (re, Co, Ci, . . . , C m ) partition of unity. 

Assume also that linear subspaces Vpj G H m (ujj) are given. These subspaces will 
be called local approximation spaces and are used to define the space 

(14) S = Sgfem ■■= { y^0i«j ; Vj e %|, 

3 

which will be called the GFEM-space. 

A basic approximation property of the GFEM-spaces is the following Theorem 
from [?]. 



6 



I BABUSKA AND V. NISTOR 



Theorem 2.2 (Babuska-Melenk). We use the notations and definitions of Defi- 
nition fK7\ and after. Let {<fij} be a (k, Cq,Ci) partition of unity. Also, let Vj £ tyj, 



y~]. 4>jVj G S, and dj = diamwj, the diameter of ujj. Then 



J 3 

\\ u - u a P \\ 2 L 2 m < kCI^\\u - Vjf^^ and 

3 



(15) C 2 \\u ~ v \\ 2 ■ 
l|V(« ~ u ap )\\l 2(n) <2 K J2( 1 , d L L2(Uj) + C ol|V(u-^)lli 2 (, 

j J 

For our method, we shall need to make some additional assumptions on the local 
approximation spaces on the covering {cjj}, and on the partition of unity {<fij}- 
We shall denote by dj the diameter of (Jj. Recall that we have assumed that ft is 
connected, the general case being immediately reduced to this one. Unless otherwise 
mentioned, we shall make the following assumptions, for some fixed values of the 
parameters A, B, Cj, n, to, <t, and h. 

It follows that 

(16) • = 1 011 Q - 

In turn, Equation l|16|l implies the second condition in Assumption A. The covering 
{ujj} will satisfy the following geometric assumption. 

Assumption A. The sets u>j are convex of diameters dj < h < 1 and there exists 
a ball u)* of diameter greater or equal ah whose closure is contained in u>j, for all 
j. Moreover, uj* C £1 if 4>j is not identically zero on Q. 

This is a non-trivial assumption, see Section El Remark 15.71 The partition of 
unity will satisfy the following condition. 

Assumption B. The family {<j>j} is a (k,Cq,Ci, . . . ,C m ) partition of unity such 
that 

(17) 4>j{x) = l iixeuj*. 

It follows from our assumptions that the sets uj* must be disjoint. Also, we have 
that 

(18) Uj\W{n) < Cihr\ l = Q,l,...,m, 

by Definition 12.11 and Assumption A. 

The local approximation spaces will satisfy: 

Assumption C. The space ^ j contains all (restrictions to u>j of) first order poly- 
nomial functions and there exists A > such that 

(19) HLff'K-) <A\\w\\ H n u . f) 
for any j, any w £ ^j, and any < I < to. 

Finally, our last assumptions is the following "inverse assumption:" 
Assumption D. There exists a constant B > 

(20) \\w\\ HHuJ]) <Bd]- s \\ W \\ HsM , 
for any j, any and any < s < t < m. 

Assumptions C and D are satisfied, for example, if j is the space of polynomials 
of degree < k ~ 1, for some fixed k>2, see Sectional 
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If A C SI is an open subset, then A satisfies Assumptions B, C, and D, but not 
A, in general. We shall hence need to single out the class of open subsets of Cl 
satisfying Assumption A. 

Definition 2.3. An open subset A C SI is called admissible if YljeJ(A) $3 = ^ 011 
A, where J{A) is the set of those indices j such that cu* C A. 

Remark 2.4. It is easy to see that A C SI is an admissible open subset if, and only 
if, A satisfies Assumption A. In particular, SI is an admissible subset of itself (see 
also Equation flUt ). Moreover, all our results on the set SI extend without change 
(including the constants) to any admissible open subset Sli C SI. 

It is proved in Section that enough admissible open sets exist if h is small 
enough. Almost all open sets below will be admissible. A typical example of an 
admissible set is obtained as follows. Fix a subset J of indices j and let G be the 
set of points where ^2j e j4>j = 1. Then the interior of G is an admissible open 
subset. 

2.2. Discrete solution. Consider the usual bilinear form 

(21) B(u,v) := / Vw Vvdx, 

Jn 

defined, for example, for u, v £ H (SI). Assume, for the purpose of this discussion, 
that CI is smooth, so that the boundary value problem Q makes sense. Then the 
solution u of Equation JTJ satisfies 

(22) B(u,v) = (g,v\ an }, 

for v smooth enough. We define then the GFEM-solution of Equation (JJJ to be 
us G S such that 

(23) B(u s ,v s ) = (g,v s \an), 

and J n us(x)dx = 0. This is possible since the form B is non-degenerate on the 
subspace So C S consisting of functions with zero integral over SI (recall that we 
are assuming, for simplicity, that CI is connected). See also Lemma 14.11 This 
also shows that we need the discretization space S to consist of functions that are 
smooth enough, which is the main reason why we are using the Generalized Finite 
Element Method in this paper. 

3. Interior estimates for the GFEM 

We continue to assume that CI is bounded and connected. Also, we do not 
assume that SI has a smooth boundary. We assume, however, that there are given 
tyj, u>* C Uj, and <pj be as in the previous section. In particular, they are assumed to 
satisfy Assumptions A, B, C, and D. This will rule out some sets SI, but any SI with 
piecewise smooth boundary enjoys this property, see Section [SJ Let S = Sqfem 
be the resulting Generalized Finite Element Space. 

All the parameters appearing in Assumptions A-D, except h, will be fixed in 
what follows (i.e., A, B, Cj, k, m, and a will be fixed). In particular, when 
we shall say that certain other constants are "independent of our choice of h and 
GFEM-space 5," we shall understand that all possible choices of a GFEM-space 
are allowed, as long as the Assumptions A-D are satisfied for an arbitrary, but 
fixed, choice of the parameters A, B, Cj, n, and a. In particular, our "constants" 
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are allowed to depend on A, B, Cj, K, to, and a. The parameter h that measures 
the degree of the refinement of our covering is allowed, however, to become as small 
as we want. Also, all constants will be independent on N, the number of open sets 
in our covering {wj}. 

All our results below remain true (with the same constants) if we replace fl by 
an arbitrary admissible open subset fii C fl. (Recall that admissible open subsets 
were introduced in Definition 12.31 ) 

3.1. if fc -approximation. We shall need a basic result on the approximation of 
functions in H k (£l) with elements in the GFEM-space S, extending Theorem 12. 21 
Only the case k = 1 will be needed in this section, but later on we shall need the 
general case. 

First, let us recall the following standard lemma. 

Lemma 3.1. Let ?pj be measurable functions defined on Q. Assume that there 
exists an integer k such that a point x G Q can belong to no more than k of the sets 
supp('0j). Let f = tpj . Then there exists a constant C > 0, depending only on 
K, such that \\f\\ 2 HHsi) < CJ2j HV'jlllrJ(n)- 

Proof. The inequality 

(24) |oi + a 2 + . • • + a M | 2 < M(\ ai \ 2 + M 2 + . . . + K/| 2 ) 

gives the desired result. □ 

Let k G Z+. We shall denote as usual 

\u\w k ^(n) '■= max ||<9 a u|| L oo (n) , ||u|| w «.,«.(n) := max ||<9 a ii|| L oo (n) , 

\a\ — k |a|<fc 

W k '°°(Q) := {u, ||tt|| w *,oo(n) < oo}, and llVa;]]^,^^ := J2j\\ d M\w<.°°(n)- In 
particular, |u|v^o,oo(q) = w a -^(n) — II u IIl°°(si)- We are ready now to prove the 
following theorem. 

Theorem 3.2. Assume that, for each j, the local approximation spaces contain 
(the restriction to u>) of the degree I — 1 polynomials. Then, for any v G H l (Q) and 
any any < k < I, there exists w G S such that 

\\v-w\\ HHn) < Ch l - k \\v\\ H , (n) 

for a constant C independent of our choice of h, S, and v G S. 

Let us notice that, by taking k — I in the above theorem, we immediate obtain 
that, using the same notation, 

(25) \H\w(a) <C\\v\\ W{n) . 

Proof. We shall use the notation and the results from [?] [Chapter 4]. Let Q y j{x) 
for any y G u>* be the Taylor polynomial of degree I — 1 at y associated to a smooth 
function / . Let then 

Wj(f)(x) := (vol^*))- 1 / Qyj(x)dy G ^ 

be the Taylor polynomial of degree I — 1 averaged over lu*. (In the terminology of 
[?], this is the "Taylor polynomial of order I averaged over lu*.) The definition of 
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Wj (/) extends to any / S H l 1 (f2) . Then we have the well known Bramble-Hilbert 
Lemma [?] [Lemma 4.3.8] 

(26) \f~^(f)U M <Ch l - k \f\ Hl(LJj)l 



with a constant C depending only on a. Let Wj — Wj(v) s 9j and w — X}j=i 0j' lt; j £ 
5. Then, using also Assumption B and Lemma 13. II we obtain 

N 

N s N s 

j — l i—1 j—l i—1 

<Cnh 2l - 2s \v\ 2 Hl{n) 

Summing over < s < k gives the desired result. □ 



Let us also record, for further use, the following well known Poincare-Friedrichs 
inequality [?], Lemma (4.3.8). (See also [?], Lemma (4.3.14), and [?], Equation 
(2.2), Theorem 14.1, and Theorem 15.3., or [?, ?].) The precise statement that we 
need is the following. 

Theorem 3.3. Using the notation of Theorem \ 2. OH we have 

(27) IMU^) < Cph\v\ H i( U:l ), 

for all v G H (u)j) satisfying f u vdx — 0, where Cp depends only on a. 
Theorems 12.21 and 13 . 31 lead to the estimate 

(28) \u - u ap \ mm < (2k)^(C 2 C p + C 2 ) 1/2 (£ l|V(u - «i)||! 2(a)j) ) 1/2 , 

provided that J (u — Vj)dx = 0. 

For k = 1, we shall need the following consequence of Theorem 13.21 which re- 
places Assumption 9.5 of [?]. Define 

(29) s< :=s GFEM nc c (n). 

That is, S < consists of the elements of Scfeai with compact support inside f2. 

Proposition 3.4. Let Qi <g Q and 9 be the distance from dfi to dVt\. Then 
there exists C > 0, independent of 9, h, and the GFEM-space S , with the following 
property. For any u £ H (Cl) with support in fii, there exists w € S < such that 

II" - wIUho) < C9~ 1 h\\u\\ H 2 {n) . 

Proof. Choose Wj and w as in the proof of Theorem 13.21 If h < 9, then Wj = 
unless ujj intersects Q,±, which gives that the closure of u>j is completely contained 
in f2. Hence the support of w constructed above is compact. For h > 9, it is enough 
to choose C large, by Theorem 13. 21 □ 
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3.2. The super-approximation property. Recall the bilinear form 

(30) B(u,v):= / Vu-Vvcte, u,v E H 1 ^). 

Jn 

Our approach follows the approach from [?] , as presented in [?] [Section 9] . See 
[?, ?, ?, ?] for related results on approximation in the "sup"-norm. Recall that 
A <s B means that the closure of A is a compact set contained in the interior of B. 
Our main goal in this section is to prove Theorem 13. 121 

Lemma 3.5. Let p be a smooth function on luj and w G tyj. Then there exists 
w £ S&j such that 

(31) \\pw - w\\m(uj 3 ) < Ch\\p\\ W 2,°o( w .)\\w\\ H i( U:j ), 

where C > may depend only on the dimension n (in particular, it is independent 
ofw, p, or j). 

Proof. We shall use the inner product induced from H (u>j). Let p € W 2 '°°(u)j) 
be given. 

To prove the lemma, we shall assume first that w € ^Sfj is a constant. Let I be 
the degree one Taylor polynomial approximation of p at the center of the ball lv*. 
Then I £ iffj, because first order polynomials are in fyj (Assumption C) and we 
have 

with C > a constant depending only on the dimension n. (This is where the 
condition h < 1 is used.) Choose w = Iw. Then 

\\pw - w\\ H i {uJj) = \\pw - Hlff 1 ^) - llP-'llwi.^^WffiO^) 

< Ch\\p\\ W 2,^ {u>]) \\w\\ H i {uJj) . 

Assume now that w S $ j is such that [w, 1) = 0, that is, w is orthogonal in 
^(ojj) to the subspace generated by constants. We then write 

P = P+P i 

where p is a constant function (say the value of p at the center of u*) and 

(32) Wp1l-M < Ch\\Vp\\ L ~ M , 

where C is a constant depending only on the dimension n. We shall choose then 
w = pw € ^ j , which makes sense since $ j is a vector space. Then 

\\pw - w\\ h i {lUj) = \\p*w\\ H i {LUj} < C||Vp*|| L = ( ^)|jw|| L 2 (a)j) 

+ C||P*IU<»K)II^IU 1 (^) ^ Ch\\Vp\\Loo( Uj )\\w\\ H i( Uj ), 

where in the last step we have used the Poincare-Friedrichs inequality for ujj (The- 
orem to estimate ||m||l2(u-) an( A Equation (|32|l above to estimate the second 
term. Here C is again a constant that may depend only on the dimension n. 
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For a general w £ 'F/, we decompose w = w\ + u>2 with w\ a constant and W2 
orthogonal to the space of constants and choose w\ and u>2 as above. Then 

\\pw -Wt- 102 1| ffi (0,^) < \\pWl ~ Wi|| ff i (a)j ) + \\pw 2 ~ W 2 \\m(u f ) 

< C , ft-||p||iy2, TC( ^ ) ||u'i|| Jf i ( ^ ) + ChWp^.co^WwuWHifa) 

< C"ft,||p|j W 2,oc (a)j) ||w|| ff l (a;j) . 

The lemma is now proved. □ 

An important technical step in our proof of the Theorem 13.121 is the following 
"super- approximation" result. 

Proposition 3.6. Let p 6 W 2,ca (Q) and w G S = Sgfem- Then there exists 
w € S such that 

(33) \\pw - 10 1| hi (n) < Cft||/9|| w a,«»( n )||u;||Hi(n), 

where C is independent of h and our choice of GFEM-space S. If u has support in 
fii <s Cl, then we can chose w to have compact support in £1 and C < C'9~ x , where 
9 is the distance from dQ to dCli and C is independent of 9, h, and S. 

As explained above, the constant C may depend, however, on the parameters A, 
B, Cj, K, m, <7, and 9, but is independent of h or of the number N of sets {ujj}. 

Proof. Let 

JV 

(34) w = ^ <f>j w j 6 & = Sgfem, Wj G 

j=i 

Let wj be the orthogonal projection of pwj onto j in the inner product of H (u>j). 
Lemma \'S . 51 then shows that 

(35) \\pwj ~Wj\\ H i {uJj) < (7/i||p||v^,=°(n)IKI|Hi(^)- 

Moreover, we have that J (pWj — Wj)dx — because the constant functions are in 
^ j and pwj — Wj is orthogonal to \& j . 

Let w :— Ylj=i ^jWj- Then \\V<j)j ||z,°°( Wj ) < C\jh by Equation ljT3)) and 

\\pwj -Wj\\ L 2 {u}j) < C P h\\pwj -Wj\\ H i( Uj) < CpCh 2 \\p\\ W 2,^ {n) \\ W: j\\ H i {u}])l 
by the Poincare-Friedrichs inequality ( Theorem 13. 3fl . and hence 

JV JV 

\\pw-wf H i {Q) = \\^<j>j(j>Wj-Wj)\\ 2 Him < (ll^illi-^^IIP^j-^illHH^-) 

JV 

+ liv^-llic.^llpwj -^Ili^) < ch 2 \\p\\ 2 w2 ^ (n) J2 HIIhm^)> 

where for the first inequality we have used also Lemma 13. II The result will follow 

now if we can prove that Y%=i \\ w j\\w(u-) - ^IMHr^n)) for an y w = ^f=i §] w ^ 
as above and C a constant independent of h and S. Indeed, we have 

JV JV 
3=1 3 3=1 
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by Assumption C (A is the constant appearing in that assumption). 

The proof of the last part is completed as in Proposition OOJ □ 

3.3. Estimates on "discrete— harmonic" functions. We shall also need the 
following "inverse property," which is somewhat similar to Assumption A. 3. in [?] 
or Assumption 9.2 in [?]. 

Lemma 3.7. There exists C > 0, independent of h, u, and the GFEM-space S 
such that 

(36) \\w\\ Hj{n) <Ch^\\w\\ 
for all < i < j < m. 

Recall that the constant m is the fixed integer appearing in Assumptions A-D. 

Proof. If i = j, we can take C = 1. Let w — J2k=i 4>k w ki with uik 6 ^k- Then 
Lemma 13. II and Assumptions A-D give 

N j 

IMlL(Q) <^^Hk\\wl,oe { u k) \\Wk\\ 2 H j-l {uk) 

k=l 1=0 

N j N 
k=l 1=0 k=l 

<C 2 h-^\\w\\ 2 L2{ny 

This proves the result for i = 0. For an arbitrary < i < j, the result follows by 
interpolation. □ 

The rest of this section follows closely the approach in the paper of Nitsche 
and Schatz [?], relying also from the survey paper [?] (which in turn is based on 
the paper by Nitsche and Schatz). There are, however, some differences in the 
assumptions that we are using, so we include complete proofs for the convenience 
of the reader. For instance, the following corollary of Lemma f3 . 71 plavs the role of 
Assumption A. 3. in the Nitsche-Schatz article [?], respectively, of the Assumption 
9.2 (Inverse assumption) in Wahlbin's article. Also, the following lemma is an 
analog of Lemma 5.2 of [?], respectively, of Lemma 9.1 of [?]. Recall that all the 
above results remain true if f2 is replaced by an admissible open subset 57i C fl. 

Corollary 3.8. Let w E S. Then there exists a constant C > such that 

\\w\\ H - Hn) <Ch^\\w\\ 

for any < i < j < m. 

Proof. If i = j we can take C = 1. Assume next that i — 0. Then 

ii ii |(w,0)| (w,w) IHLl=(o) x in ., 

\\w Ur-i(n) = sup-|r-|i > = > C h J \\w\\ L 2 {n) , 

\m\Hi(n) \m\HJ(n) \\w\\hhu) 

by Lemma 13.71 For the other values of j, the result follows by interpolation. □ 

We shall denote A <s B if A, the closure of A in R 2 , is a compact subset of the 
interior of B. Also, we shall denote 

(37) S < (A) := {u G Sgfem, supp(u) € A} 

for any admissible open subset A C fi. In particular, S < (fl) = S < . 
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Lemma 3.9. Let A <e A\ d be admissible open sets and 9 — dist(dA,dA\). 
Then there exists C > 0, independent ofh, S, A, and A\ with the following property. 
IfweS — Sqfem and 

(38) B(w,x)=0, forall X eS < (A 1 ), 

then, for h small enough, \\w\\jjir A -\ < CHidH^^), with C depending only on 9 
and not on A, A\, h, or S. 

Proof. The proof is the same as the one in [?, ?], using Lemma 13.71 in place of 
Assumption A. 3, respectively Assumption 9.2 ("Inverse assumption"), and Propo- 
sition l3~^l in place of Assumption A. 2, respectively Assumption 9.1 ( "Superapprox- 
imation"), of [?], respectively [?]. The constants "C" below are allowed to depend 
on 9. 

Let us chose A <g Aq <<= A\ admissible open sets such that the distances between 
the boundaries of these sets are > 9/C. This is possible if h < c9 (see Section EJ. 
Also, let uj e C c °°(Ao), with u = loni,w>0on A , and ||a;||virt.~(n) < C9~ k , for 
k = 0,1. Then, by Equation (|38Jl . we obtain 

\Nw\\ 2 L 2 (A) < (Vm,wVw) = (Vw,V(ww)) - (Vw,V(w)w) 

= (Vw, \7(uiw - %/})) + i(w, (Au)w), 

where the inner products are in L 2 (A ) and tJj G S < (A ). Proposition 131)1 then 
gives 

||Vw||£ 2(A) < Ch\\iv\\ W 2,~ { n)\\w\\ 2 HHAo) + C\\w\\I HAq) , 
which, in turn, implies 

(39) \\w\\m { A) <C(ft 1/2 ||HI H 1 (A ) + \\w\\l2{A ))- 

We now repeat the argument for Aq <e A\ <s 51 (and A replaced by Aq and Ao 
replaced by A\), which gives 

(40) \\w\\ HHAo) <C(hV 2 \\w\\ 

Combining Equations l)39|) and (|4U[) and using also h < 1, we obtain 

(41) \\M\m(A) < C(h\\w\\ H i {Al) + C\\w\\ L 2 (Al) ). 

Since A\ is admissible, we can use Lemma |3.7I with SI replaced with Ai to obtain 
^IMI-ff^Ai) < C|l u; lli 2 (Ai): with C independent of h, w, and the GFEM space 5* 
(as long as S satisfies the Assumptions A-D). Then 

(42) \H\m { A)<C\\w\\ L 2 (Al) . 

The proof is now complete. □ 

We shall need the following simple estimate. 

Lemma 3.10. Let §(x) = log | a; if n = 2, $(x) = |x| 2_ ™ 7 if n ^ 2. Let U be a 
fixed bounded open subset o/R™. Then there exists C > 0, which depends only on 
U , such that 



$ * u{x) := / $(x — y)u(y)dy 
Jy 

satisfies 

||$ * u\\ H i+2 {U) < C\\u\\ H i {u) , 
for any I G R and any u G C£°([/). 
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Proof. Let ui e C C °°(R") be equal to 1 on U. Then Tu{x) = w(x)[$ * (wu)](x) is 
a pseudodiffcrential operator of order —2 with compactly supported distribution 
kernel. Hence it is bounded as a map H l (M. n ) — > H l+2 (W l ) [?, ?]. The statement 
follows by restricting T to C C °°(C/) C H l (R n ). □ 

We define 

(43) || U || =sup I^L< N | O^veC^(A) 

° \\ V \\H'(U) 

for any open set U, any u G L 2 (U), and any s > 0. We define Hg S (U) to be the 
completion of L 2 (U) in the norm ||u||#~-«/m- Then Hq S (U), s > 0, identifies with 
the dual of # S (C/). 

We now prove the following lemma. 

Lemma 3.11. We keep the notation and assumption of Lemma \3.!A In particular, 
we assume that w 6 S satisfies Equation JSSJ. Then, for h small enough, 

(44) H\ LHA) <C\\w\\ H - m{Al) , 

where C is a constant depending only on the distance 9 from dA to dAi and not 
on A, A\, h, or S. 

Combining Lemmata 13 . 91 and 13 . 1 II we obtain 

(45) \H\hhA) < C\\w\\ H - m(Al) , 

for h small enough and any w € S satisfying the assumptions of Lemma l3.9l 

Proof. Let A <s B (g P>i <g A\ be such that the distances between the boundaries 
of these sets are > 9/C. This is possible if ft, < c9. Note that by Lemma f3. 91 we 
have 

(46) IMIffH*) ^ C W w \\l-{a,). 

For any v 6 C£°(A), let V := c„$ *«£ H l+2 (Bi), where c„ is chosen such that 
AV = v (see [?]). Lemma f3. 101 then gives 

(47) \\V\\ H i + ! iBl ) < C\\v\\ H i {Bl) = C\\v\\ H i {A) , l£Z + , 

for some constant C that depends only on B\. Let oj S C^°(Bq) with w = 1 on 
A and ||w||w r »t.»(n) < C9~ k , for < k < m. Since <E C%°(B ), we know from 
Proposition 13 . 41 that there exists x G S < (Bi) such that 

(48) - xIIh^Bi) ^ C/i||o;^||^ (Bl) < Ch\\V\\ m{Bl) < Ch\\v\\ L , [A) . 
Then, for any v 6 C c °°(^), 



(tOji;)^ = (ww,v) A 



/ wwAydx = / ujwAVdx = / V(ww) • WVdx 

J A JB JBq 

-■ / w(2Vuj-VV ~VAuj)dx+ / Vw • V(u/F - 

«/ Bq J B\ 



for any \ G S < (B\), where the inner products are calculated on the indicated 
sets. Then, by combining Equations (|4T)1) . ljT7l) . and as well as Lemma 1X^1 we 
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obtain, for all / > 0, 

\(w,v)a\ < C\\w\\ h -i-i {Bq) \\V\\ h i+2 {Bi) +/i||w|| J? 1( Bi )II^IIh 2 (B 1 ) 

^ C (W W \\h- 1 - 1 (A 1 ) + h \\ w \\Li(A 1 ))\\v\\ H '(A), 

and hence 

(49) \\ w \\h- 1 (A) ^ CtHHIff-'- 1 ^!) + h \\ w h 2 (A 1 ))- 

Let us now choose a sequence of open sets A (<= B\ (e B 2 <s . . . (s B m d fl with 
all distances between the boundaries greater or equal c8. Changing notation and 
iterating Equation 1)49(1 . we obtain 

(50) \\w\\ L 2 {A) < C(\\w\\ H -i (Bi) + h\\w\\ L 2 [Bl) ) 

< C(\\W\\ H -2 {B2) + h\\ W \\ L 2 {B2) ) <... 

^ C (W W \\H-™(B m ) + h \\ w \\^(B m ))- 

We now repeat the above reasoning. We change notation again, so that, this 
time, B m becomes B\, then we chose as before a sequence of open sets 

A € Bx <s B 2 <s . . . € B rn <e O 

with all distances between the boundaries greater or equal cO. Then we iterate 
Equation ((50(1 . and obtain, 

IMk 2 (A) < C(\\w\\ H - m{Bi) + h\\w\\ L 2 {Bl) ) 

< C(IHI* -« (Ba ) + ^IMIl 2 ( b 2 )) < ••• < C([[«;|| ff -« (Bm) + ftlHU'dw) 

^ C (ll U, lli?-"'(B m ) + \\ w \\H-™(B m )) < C\\w\\ H - m{Bm) , 

where at the end we have used the inverse property ||u>||.z,2(£/) < h~ m \\ w \\ H- m (u) f° r 
any admissible open set U <Z ft (see Corollary 13.8(1 . The proof is now complete. □ 

3.4. The interior error estimate. The following result, the main result of this 
section, is an analog of [?] [Theorem 5.1] and of [?] [Theorem 9.2]. 

Theorem 3.12. Let A <s= B C Q be admissible open sets. Then there exists C > 
with the following property. IfuG H 1 ^) and us G S are such that B(u—us,x) = 
for all \ e S < := S < (£1), then for h small enough, 




The constant C depends only on the distance 8 = dist(dA, OB) and not on h, S , 
or the sets A and B . 

Proof. Let A <s A x g A 2 <s B c ft. Choose ui G C^°(A 2 ) such that ui = 1 on Ax. 
Let Px be the iP(f2) orthogonal projection onto S < (A X ) C S C i^ 1 (Sl). Then on 
Ai 

(51) u — us = {<jJu — Pi (wu)) + (Pi (uu) — us) • 
Then, by the general properties of orthogonal projections, we have 

(52) \\wu - Px(uu)\\ H i (n) < \\uju\\ H i {n) < C\\u\\ H i {B} . 
Hence 

(53) \\uju - Px(u}u)\\ H -™(Ai) < - Pi(uu)\\m(n) < C\\u\\ H i( B y 



16 



I BABUSKA AND V. NISTOR 



Let w = P\(u>u) — us- Then B(w, \) = BIluu — us, x) — B(u — us, x) = 0> f° r 
all x G S' < (j4i ) , and hence w satisfies the assumptions of Lemmata 13.91 and 13. 1 II 
From this, using also Equations 1)45(1 . (|5ip. and (|53|l . we obtain 

( 54 ) IMIff 1 ^) - C 'll w llH- m (A 1 ) < - Pi{uju)\\ H -m {Al) + \\u - u s \\ H -™(Ai) 

< \\u\\m(B) + \\u-u s \\h-™ (AO- 
Equations H51H54JI then give 

\\u-us\\ H - m (A) <\\uu- P!{uju)\\ H - m(A} + \\w\\ H - m{A) by (EU 

< C||u|| H i(B) + II" - u s \\ H - m(Al) by dS3-(|HJ- 

The desired result follows by replacing u and us with u — x and, respectively, us~x, 
with x in S = Sqfem- D 

4. Approximate solution of the Laplace equation with distribution 
boundary conditions using the gfem 

We shall consider the same setting as in the previous sections, and we shall 
further assume that fl is a smooth domain. In particular, S will be the GFEM- 
space associated to a partition of unity {4>j} subordinated to the covering {loj} 
of f2 and local approximation spaces {^j}- We shall continue to assume that 
Assumptions A-D are satisfied, for a fixed choice of constants A, B, Cj, k, m, and 
a. In particular, the "constants" below are allowed to depend on these parameters, 
but are not allowed to depend on h or the specific choice of the GFEM-space S, as 
long as the constants above remain the same. 

We shall denote by v the outer unit normal vector to <9f2. By d v u(x) we shall 
denote the directional derivative of a function u in the direction of v, at some point 
x on the boundary. 

Let u € H 1 ~ k (fi). We want to make precise in what sense we shall say that 
"Au = as a distribution on fl" Recall first that the space H 1 ~ k (Vl), k € Z + , 
k > 1, was defined as the dual of H k ~ 1 (U) (see Section^]). We shall write (u, v) := 
u{v) e C for any u £ H 1 ^ k (fl) and v £ H k ~ 1 (tt) for the "value of u evaluated at 
v." We can hence define by duality dj := -d* : H 1 ~ k (Q) — ► H~ k (ft). This leads 
to a definition of Au £ H~ 1 ~ k (Q), for any u £ H 1 ~ k (fl) by 

(Aw, v) := (u, Av), for any v £ H 1 ~ k (Q). 

However, this turns out to be too strong a condition. Instead, we shall require 

(55) (Au, <j)) := (u, A<j>) = 0, for any £ C c °°(fi). 

We shall say that Am = as a distribution on fl whenever Equation (|55|l is 
satisfied. 

Typically, u as above will arise as a solution of a boundary value problem, for 
example, as a solution of the boundary value problem lfT|l. In [?,?,?] it was 
explained how to define the traces (or restrictions) u\en and 8 v u\qq, for any u £ 
H 1 - k (fl) satisfying Au = as a distribution on Q. More generally, we define 
u\oa £ H l l 2 ~ k {VL) and d u u\on £ H~ 1 l 2 ~ k (Q) by linearity, for u = u\ + u 2 , where 
Mi £ H 1 ~ k (£l), Ami = as a distribution on fl, and M2 £ i/ 1 (f2). We use this to 
define B(u, v) by 

(56) B{u,v):=-(u,Av) + (u\ dn ,d„v), for any v £ H 1+k (n), 



DISTRIBUTIONS 



17 



for u = ui + U2 as above. In view of Green's formula (see [?], for example), 
B(u,v) = Jo Vu • Vvdx if u,v € H 1 ^), as originally denned. It is not clear how 
to define -B(m, for arbitrary zt G H 1 ~ k (£l), since (tt, Aw) is denned but the traces 
of u may not defined in general. 

From now on, we shall fix u such that, and 

(57) Au = 0, ueH 1 -*^), 

as a distribution on fl, where k G Z+, m + 1 > k > 0, is also fixed. We do not 
assume that Ait = in i? _1_ ' c (ri) (i.e., we do not assume (u, Av) = for all 
v G we only assume (u,Av) = for all £ C~(fi)). We shall also 

assume that 

(58) (u, 1) = (d v u\ e n, 1) = 0. 
We have the following. 

Lemma 4.1. There exists a unique us G S such that (us, 1) = and 

(59) B(u-u s ,v s ) = 0, 
for all vs G S. 

Proof. Let So be the subspace of the GFEM-space S consisting of functions %o G S 
with J^x^da; = 0- The bilinear form B is non-degenerate on So- This gives the 
existence of a unique us G So such that Equation i|59[) is satisfied for all vs G So- 
Since S = S + C, the result follows from B{u, 1) = B(u s , 1) = 0. □ 

We also have the following simple estimate. 

Lemma 4.2. With u as in Equation Ij57(l above, we have 

\B{u,v)\ < C\\u\\ H i-k {n) \\v\\ H i+k {n) , 

for any v G H 1 ~ k (£l) and a constant C depending only on il. 

Proof. By definition, we have 

(60) \B(u,v)\ = \-(u,Av) + {u\ d n,d v v)\ < \(u, Av)\ + \{u\ d n,d„v}\ 

< \H\m-"(Q)\\^\\ H -l+k^ + \\u\\ H l/2- k{dn) \\v\\ H -l/2 + k {m) 

< C\\u\\ H i- k{n) \\v\\ H i+ k{n) . 

This completes the proof. □ 

We continue with more lemmata. Recall that k G Z+. 

Lemma 4.3. We have \\us\\h 1 {si) ^ Ch~ *\\u\\gi-k(Q\ for a constant C depending 
only on f2. 

Proof. The Poincare-Friedrichs inequality and Lemma 14 . 21 give 
\\us\\ 2 H i(n) < CB(u s ,u s ) = CB(u,u s ) < C\\u\\ H i-h^\\u s \\ H i+h^ 

< Ch~ k \\u\\ H i- k ^\\us\\m(n), 
where in the last inequality we have used Lemma 13.71 □ 

This gives the following corollaries. 

Corollary 4.4. We have \\us\\H 1 - k (fi) < ^IMIfl^-fcfn) for a constant C depending 
only on fl. In particular, \\u — ws||iji-fe(n) < C||' u ]].H" 1 -'«(n) • 
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Proof. The result is well known for k — since us is the i?-orthogonal projection 
of u onto S (see also Cea's Lemma, [?, ?]). We shall therefore assume that k > 1. 

Let v <E H k ~ 1 (£l) be arbitrary. Let c e C be such that J n (v — c)dx = 0. Then 
we can find V 6 H k+1 (fl) such that 

-AV = v-c, [ Vdx = 0,d u V = 0, and \\V\\ Hk+ i m < C\\v\\ Hk -i {n} , 
Jn 

where C is a constant depending only on Q. Also, chose w £ S such that 

\\w\\H"+Hn) < C\\V\\ H k+i {n) and ||U - Hlff^n) < Ch k \\V\\ H k+i {n) . 
This is possible by Theorem 13. 21 Then 

(u s ,v) = (u s ,v- c) = ~(u s ,AV) = -{u s ,AV) + {u s \ao,d v V) = B(u s , V) 
= B(us, w) + B(us, V — w) — B(u, w) + B(us, V — w). 
Using also Lemmata 14.21 and 14.31 this gives 

\{us,v}\ < C\\u\\ H i-k {n) \\w\\ H i+k {n) + \\u s \\ H i m \\V - w\\ H i {n) 

< C\\u\\ H i- k{n) \\V\\ H k + i (n) +Ch^ k \\u\\ H i- k{n) h k \\V\\ H k+i {n) 

< C[[«|| H i-*( n ) 11^11^+1(0) < C\\u\\ H i- k{n) \\v\\ Hk -i in) . 

This gives the result since \\u s \\ H i-k^ n) := sup \(u s , v)\/\\v\\ H k-i^, v 0. □ 
Similarly, 

Corollary 4.5. We have \\us\an \\H 1 / 2 - k (dn) < C|l u ll_H' 1 - fc (o) f or a constant C 
depending only on il. In particular, \\(u — us)\dn\\H 1 / 2 - k (an) ^ C|l u ll_ff 1 - fc (n) • 

Proof. The proof is similar to that of the previous corollary. Let v 6 H~ 1 / 2+k (dVl) 
be arbitrary. Let c £ C be a constant such that J af2 vdS = J n cdx. Then we can 
find a unique W £ H 1+k (£l) satisfying 

AW = c, f Wdx = Q, d v W = v, and \\W\\ HW{n) < C\\v\\ H -i /2+Hdn) , 
Jn 

for a constant C > depending only on fl. Using also Theorem l3.2l we choose w E S 
such that |H|ff*.+i(n) < C||PU|| ff fc+i ( n) and \\W - w\\ H i {n) < Ch k \\W\\ H k+i {n) . 
Then, using also {us, 1) = 0, we obtain 

(«s|en, v) = (uslan, = («s, AFU) + W) = B(u s , W) 

= B{us, w) + B(us, W — w) — B(u, w) + B(u s , W — w). 

Using Lemmata 14.21 and 14.31 we then obtain 

\{us\dn,v}\ < C\\u\\ H i- k{n) \\w\\ H i+ k{n) + \\us\\ H i{a)\\W - w\\ H i {n) 

< C\\u\\ H i- h(Q) \\W\\ H x+k m +Ch- k \\u\\ Hl - k{n) h k \\W\\ H i+k m 

< C\\u\\ H i-k(n)\\W\\Hi+k(Q) < C||«||H l -*(n)[|^||fl--V2+*(8n)) 
which completes the proof in view of the definition of ||its||iji/2-fcm). □ 

We complete our sequence of estimates with the following result. 
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Proposition 4.6. Let k + 7 < m + 1, k, 7 G Z+. Assume that each local approxi- 
mation space j contains the polynomials of degree l + fc + 7. Then the error u — us 
satisfies 

\W - Us||ifi-*-r(fi) < C7l 7 ||u|| ff i- fc(n) , 

ura'i/i a constant C independent of h or the GFEM-space S , but possibly depending 
on f2 and the parameters A, B, Cj, n, a, and m. 

Proof. Let v G H~ 1+k+ ~ < (£1) be arbitrary. Let c be a constant such that (v — c, 1} = 
0. Then there exists a unique V G H 1+k+ ~ f (fl) such that 

-AV = v-c, [ Vdx = Q, d„V = Q, and \\V\\ H i +k+Hn) < C\\v\\ H - 1+k+Hn) , 
Ja 

for a constant C > depending only on f2. 
Then, for any w *E S, 

(61) (it — us, v) — (u — us, v — c) = —(u — us, AV) 

= -{u- u s , AV) + ((« - u s )\gn, d v V) = B(u - i* s , V) = B(u - u s , V - to) 

= -(it - us, A(V - to)) + ((it - us) I du, d v w). 

Using Theorem l3.2l we chose w G S such that ||u>||#i+k+ 7 (Q) < C|| V|| ^l+fc+T (n) an d 
||V-tu||jji+k(n) < Ch~t\\V\\ H i+k+-, (n y In particular, 

\\d v w\\ H -i/2+k {dn) = \\d v (V - w)\\ H -i,2 +k(an) < h 1 \\V\\ H i+k+~ /in) . 

From ||V||ifi+fc+-r(n) < C||u|| ff -i+fc+ T (n), Corollaries 14 . 41 and 14.51 and Equation (|ol"Jl . 
we then obtain, 

\(u-us,v)\ < ||u-us|| ff i-fc(o)||U-u;|| H i +fc(a) + ||u-us|| ff i-fc(o)||^w|| H -i/2 + )= (an) 

< C/i 7 ||u|| H i- fe(fi) ||u|| ff -i + (=+^ ( n)- 
The proof is complete. □ 

Theorem 13.121 and Proposition 14.61 then give the following result, which is the 
main result of this paper. 

Recall, for the following theorem, that S = Sqfem is the GFEM-space associ- 
ated to a partition of unity satisfying Assumptions A-D. Also, recall that we have 
fixed u G H 1 ^^) satisfying (u,A<t>) = for all cj) G C™(fl) and that us is the 
GFEM-approximation of u (i.e., given by Lemma l4.1|) . 

Proposition 4.7. Assume the local approximation spaces tyj contain the polyno- 
mials of degree k + 7 + 1 and let A <s B <g Q be admissible open subsets. Then for 
any —l + k< —l + k + "f<m and any I > 1, k, 7 G Z+, we have 

\\u - usWh^a) < Ch l \\u\\ H i+i( B) + C7i 7 ||u - us\\h-i+*+~i(b)- 

Proof. This follows from Theorem 13.121 and from 

in| II" - x\\m(B) < Ch l \\u\\ H i+i {B} . 

□ 

By taking I — 7 and using also Proposition 14.61 and Equation (J2J), we obtain 
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Theorem 4.8. Assume the local approximation spaces \& j contain the polynomials 
of degree k + 7 + 1 and let Aq <s f2 be an admissible open subset. Then for any 
— 1 + fc < — 1 + fc + 7 < to, fcj £ Z +; we have 

\\u-usWhhao) < C/i 7 ||u|| H i-fc ( n)- 
The constant C above is independent of h and S , but may depend on the parameters 
A, B, Cj, k, a, I, 7, and to, as well as on 9, the distance between dA a and dfl. 

5. Polynomial local approximation spaces 

In this section we shall verify that the Assumptions A-D are verified if we choose 
to be the space of (restrictions to Wj of) polynomials of degree < m, m > 1. 
Most results of this section are either elementary or well known. We include them 
nevertheless for the benefit of the reader and for completeness. 

In this section, the set of polynomials of degree to will be denoted Q m . Also, for 
any ball B of radius r, we shall denote by tB the ball with the same center as B 
and radius tr. 

Lemma 5.1. There exists a constant C > 0, depending only on n, to, and t > 0, 
such that for any ball B C M" and any Q £ Q m , we have \\Q\\L 2 (tB) < C||QIIl 2 (b) • 

Proof. For any fixed B, Q 1— » ||Q||L 2 (ts) an d P 1— > ||Q||l 2 (b) are two norms on 
the finite dimensional space Q m of polynomials of degree < m, and hence they 
are equivalent. This gives the result, except the independence of C on B. But all 
balls are affine equivalent and the L 2 -norm is scaled by the (square root of the) 
determinant of the matrix of the affine transformation. Thus the constant C can 
be chosen to be the same for all balls B. □ 

This gives immediately the following corollary. 

Corollary 5.2. There exists a constant C > 0, depending only on t > 0, n, 
and m, such that for any ball B C M™ and any polynomial Q £ Q m , we have 
\Q\w(tB) < C\Q\ H i (B) and \\Q\\ H '(tB) < C\\Q\\ H i (B) , 0<l <m. 

Proof. Use Lemma f5. II for all derivatives Q^ a \ where a is a multi- index such that 
H = l or \a\ < I. □ 

We now establish to the following "inverse property." 

Lemma 5.3. There exists a constant C > 0, depending only on m, a, and n such 
that ||Q( q )||l 2 (B) < Cr /- I Q I \\Q\\h 1 (b) f or an V I < \ a \ < m , an U Q 6 Qm, and any 
ball B of radius r. 

Proof. Let us prove first the result for I = 0. That is, we need to prove that 
\\Q {a) \\LHB)<Cr-\ a \\\Q\\ L2{B) . 

Let B\ = -Bi(O) be the unit ball centered at 0. Then Q 1— > HQ*-"-* ||l 2 (Bi) is a 
semi- norm on Q m , the space of polynomials of degree at most to, and hence it is 
bounded by the norm Q ^ \\Q\\ L 2 {Bl) . Thus \\Q^ \\ L 2 [Bl) < dWQW^^y Let L 
be an affine transformation mapping B\ onto the ball B of radius r consisting of 
the composition of a translation and a dilation of ratio r. Then 

\\Q^\\ LHB) = det{Lfl 2 \\Q^ o L\\ LHBl) = det{L)^ 2 r-^\\(Q o L)^\\ mBl) 
^ddet^^r-WHQoLlU^) = Cr-H||Q|| L2(B) , 



DISTRIBUTIONS 



21 



for any Q e Q m . 

Assume now that \a\ > I > 0. Choose (3 < a, |/3| = I. Then 

Wrr-PDPQWvW < Cr-\ a -^\\D^Q\\ L2(B) < Gr , -I°'I||Q|| H , (B) . 

This completes the proof. □ 

The relevant "inverse property" now follows. 

Proposition 5.4. There exists a constant C > 0, depending only on a, m, a, and 
n such that \\Q^ |U 2 (n) < CV'~I q| for any I < \a\ < m, any Q e Q m , 

any ball B of radius r , and any convex set fl contained in o~ l B . 

Proof. This follows from Corollary 15 . 21 and Lemma f5. 31 □ 

We now prove the following elementary lemma. 

Lemma 5.5. For any A C= f2 and any k 6 Z+, we can construct admissible open 
sets A =: B <s B x <e B 2 <& . . . <e B k <e B k+1 := fl such that Cdist(dB j ,dB : j) > 
9/k, where 6 :— dist(dA, d£l) , provided that Ch < 6, where C depending on n only 
(in particular, C is independent ofh). 

Proof. Take k = 1, for simplicity. The general result is proved similarly or by 
iterating this case. Let C — 4. Let U be the union of all open sets cuj at distance 
at most 9/ A from A. Let J be the set of indices j such that <f)j ^ on U and let G 
be the set where Y^je J = ^" ^ e then define B\ to be the interior of G. □ 

5.1. Partition of Unity. We show in this section that, for a suitable set f2, we can 
choose a family of partitions of unity {4>j}, subordinated to the covering {tOj}, with 
h — > but with all the other constants fixed. The proof of the following theorem 
is not constructive. A constructive proof, suitable for numerical implementation, 
will be included in a forthcoming paper where we will also discuss the numerical 
implementation of the GFEM for boundary value problems with distributional data. 

Theorem 5.6. Let Q be a bounded open set with piecewise C 1 -boundary. Then 
there exist constants A, B, Cj, k, m, and a, such that, for any small enough h > 0, 
we can construct a partition of unity {<fij} subordinated to the covering {ujj} and 
satisfying the Assumptions A-D of Section (J2J for the given value of the parameters 
A, B, Cj, k, m, a, and h. 

Proof. Let us first construct the covering {uij} and the subsets {w*}. The index 
j will belong to a set of points of fi, the centers of those balls. (So uij and uij will 
have the same center, namely j.) 

Let e > be small enough such that, for any y 6 satisfying dist(y, dil) < e, 
there exists a unique z € dQ with dist(y,dtt) = dist(y,z). Let r r C f2 be the 
set of points at distance r, r < e/2 to <9f2. Then T r will be a piecewise C 1 curve, 
bounding a domain diffeomorphic to f2. Choose on T r a maximal set of points Xh 
containing the vertices of T r and at distance at least r/2 from one another. The 
maximality of then guarantees that the distance between any two consequtive 
points in Xh is at most r. Let 

f(r) := sup dist(y,T r ). 
yedn 

A geometric argument based on the assumption that dft is C 1 then shows that 
lim/(r)/r = fi, as r — > 0, with [i finite. Let then a^ 1 > max{4(^ + 1),8} and 



22 



I BABUSKA AND V. NISTOR 



h := ra^ 1 /2. This guarantees that the balls with centers in Xh and diameter h will 
cover the region between L and dft. Let Yh be a maximal subset of f2 containing 
Xh and such that any point in Yh is at distance at least r to the boundary dVt and 
at distance at least r/2 from any other point in Yh. The choice of Yh shows that the 
balls io* of diameter ah = r/2 and center at the points Yh will be disjoint, whereas 
the balls u)j with the same centers and diameter h will cover the interior of T. 

Let us now construct the partition of unity (fij. Let 1 > ip(t) > 0, t > be a 
smooth function such that ip(t) = 1 if t < a and ip(t) = if t > 1. Also, let > 
be a smooth function such that = if t < 1 and 1 > > if t > 1 and 
C(t) = 1 if t > 2. Define then 

(62) Vj (x) - ip(dist(x,j)/h) J| ((dist(x,j')/h). 

Let f?(x) = arL( i = 7 ?i( a; )/ ? ?( a ')- We observe that the number of 

factors 7^ 1 in the definition of i]j is bounded by a constant independent of x, j, 
and any of the choices above. Therefore {0 j} is our desired partition of unity. □ 

Some assumptions on the domain in the above theorem are necessary, as shown 
by the following remark. 

Remark 5.7. The non-Lipschitz domain 

fl c := {(x, y), -x 2 < y < x 2 , x 2 + y 2 < 1, x > 0} 
will have no covering {uij} satisfying the Assumptions A-D. 
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